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What is claimed is: 

1. A method of using a processor to analyze data signals representing 
tomography scan images of an organic object, the method comprising: 

a) receiving data representing successive tomography scan 
5 images of said object; 

b) performing calculations by imposing an inequality constraint to 
determine dynamic data values from said data, each of said 
dynamic data values representing a physical property of said 
object at a respective corresponding one of a plurality of voxels 

10 of said object at a respective corresponding time; and 

c) producing a representation of said dynamic data values, 
representing said physical property at said voxels at said times. 

2. A method as claimed in claim 1 wherein performing said calculations 
includes minimizing a figure of merit function subject to said inequality 

15 constraint. 

3. A method as claimed in claim 2 wherein performing said calculations 
includes minimizing said figure of merit function subject to linear basis 
functions constrained by a linear inequality constraint selected from the 
following group wherein N represents the number of said respective 

20 times: 

a) X., > x,^ > x,3 . . . > > x,^, > 0 ; 

b) 0 < X,., < x^^ <x.^...< < x.^ ; 

c) 0 < X,, < - • < x.^_^ < x,^ > x,^^^ . . . > x,^_j > x,^ > 0 ; 

d) 0 < x,.j < x,.2 • • • < X. and x. . . . > x,.^_, > x. ^ > 0 ; 
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e) 0 > - ^ -Xj^ >...> Xi^^_i - jjc,.^ and 
X,, >0,...Xij, >0; 

f) 0 < - Xi2 ^ - -^(3 ^■■■^ - X;^^ and 
Xn >0,...X;^ >0; and 

g) 0 > - ^ x,.2 - x,.3 > . . . > x,.^ - X;^^j and 

A method as claimed in claim 1 wherein receiving data includes 
receiving data representing successive images detected at a Single 
Photon Emission Computed Tomography (SPECT) imaging device. 

A computer-implemented method of producing a representation of a 
measurable property which varies in time and space, the method 
comprising: 

a) receiving a plurality of sets of values representing 
measurements of said property across an object at respective 
measurement times, each set being associated with a 
respective measurement time; and 

b) producing a plurality of sets of values representing said property 
at a plurality of locations throughout said object at said 
respective measurement times, by minimizing a figure of merit 
function relating said values representing said measurements 
with said values representing said property at said plurality of 
locations, with a shape constraint Imposed on said values 
representing said property at said plurality of locations. 

The method claimed In claim 5, wherein minimizing said figure of merit 
function further comprises solving for numerical values describing said 
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property for each of said locations throughout said object, said 
numerical values being constrained by said shape constraint. 

7. The method claimed in claim 5, wherein minimizing said figure of merit 
function further comprises solving for linear basis functions describing 

5 said property for each of said locations throughout said object, said 

linear basis functions being constrained by said shape constraint. 

8. The method claimed in claim 5, wherein minimizing said figure of merit 
function further comprises minimizing a sum of squares of a difference 
between a product of a linear operator and said values representing 

10 said property, and said values representing said measurements. 

9. The method claimed in claim 8, wherein minimizing said figure of merit 
further comprises solving the following equation: 

a) 

b) 

15 

c) 



20 d) 



i) X,, >x,, >x,3...>x.^_j>x,^>0; 



Minimize Q^a jk~'{llCy,x,(h)-^jk^ 



, wherein: 



djk comprises said plurality of sets of values representing 
said measurements dj of said property across an object at 
respective measurement times tk, 

Xi{tk) comprises said plurality of sets of values representing 
said property x at a plurality of locations i throughout said 
object at said respective measurement times tk, 

x^tj) = Xik are numerical values subject to said shape 
constraint comprising a linear inequality constraint selected 
from the following set of constraints, wherein N represents 
the number of said measurement times: 



25 



ii) 0<x,, <x,2<x,3...<x,^_^<x,^; 
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iii) 0 < X,, < X,, • ■ • < < x,^ > x,^^^ . . . > x.^_, > x.^ > 0 

iv) 0 < x,.j < • • • < _^ and 
Xj._^^^ . . . > x^- > x^^^ > 0 ; 

V) 0 > X., - X,-2 > X^2 - ^ • • ■ ^ - 

x,i>0,...x,^>0; 

vi) 0 < X;j - X;2 < x,.2 - x,.3 < . . . < x,.^_j - x,.^ and 
x,-! > 0,...x.^ > 0; and 

vii) 0 > x,.j - X;2 ^ x,.2 - x,.3 > . . . > x,-^ - x,-^^i and 

e) Cp and a^jk are weighting factors; and 

f) R(x) is a regularizing temn. 

The method claimed in claim 5, wherein minimizing a figure of 
merit function includes executing a math program on said values 
representing said measurements, said shape constraint being pre- 
specified to said math program. 

The method claimed in claim 5, further comprising producing a 
graphical representation of said object, for controlling a display device. 

The method claimed in claim 11, further comprising producing a time 
varying graphical representation of said object to represent a change of 
said property over time. 

The method claimed in claim 1 1 , further comprising transmitting said 
graphical representation to a display device for displaying a 
representation of said property. 
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14. The method claimed in claim 13, further comprising displaying said 
representation of said property. 

15. The method claimed in claim 14, further comprising producing 
successive images representing said property at successive instants in 

5 time, in response to said representation of said property. 

16. The method claimed in claim 16, further comprising producing said 
successive images to depict a three dimensional representation of said 
property in said object which varies according to changes in said 
property over time. 

10 17. The method claimed in claim 5, further comprising producing said sets 

of values representing measurements of said property across said 
object at respective measurement times. 

18. The method claimed in claim 17, wherein producing comprises 
measuring radioactivity across said object. 

15 19. The method claimed in claim 18, wherein measuring radioactivity 

includes operating a Single Photon Emission Computed Tomography 
(SPECT) imaging device. 

20. A signal produced by the method of claim 5, said signal representing 
said plurality of sets of values representing said property at a plurality 

20 of locations throughout said object. 

21 . An apparatus for producing a representation of a measurable property 
which varies in time and space, the apparatus comprising: 

a) a receiver for receiving a plurality of sets of values representing 
measurements of said property across said object at respective 
25 measurement times, each set being associated with a 

respective measurement time; and 
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b) a processor circuit in communication with said receiver, said 
processor circuit being programmed to produce a plurality of 
sets of values representing said property at a plurality of 
locations throughout said object at said respective measurement 
times, by minimizing a figure of merit function relating said 
values representing said measurements with said values 
representing said property at said plurality of locations, with a 
shape constraint imposed on said values representing said 
property at said plurality of locations. 

The apparatus claimed in claim 21, wherein said processor circuit is 
programmed to solve for linear basis functions describing said property 
for each of said locations throughout said object, said linear basis 
functions being constrained by said shape constraint. 

The apparatus claimed in claim 21, wherein said processor circuit is 
programmed to minimize a sum of squares of a difference between a 
product of a linear operator and said values representing said property, 
and said values representing measurements of said property. 

The apparatus claimed in claim 23, wherein said processor circuit is 
programmed to solve the following equation: 



a) mn\m\ze }^^(j.~H^Cy,x,(t,)-dj,] wherein: 



b) djk comprises said plurality of sets of values representing said 
measurements dj of said property across an object at 
respective measurement times t/c', 

c) Xi(ti() comprises said plurality of sets of values representing 
said property jc at a plurality of locations i throughout said 
object at said respective measurement times tk. 
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d) Xi(tO =Xik are numerical values subject to said shape constraint 
comprising a linear inequality constraint selected from the 
following set of constraints, wherein N represents the number 
of said measurement times: 

i) X,, >x., >x,3...>x,.^_, >x,^ >0; 

ii) 0 < < < . . . < < x-f^ ; 

iii) 0 < < • • • < < x,^ > x,^^^^ . . . > X,. > x,^ > 0 

iv) 0 < x,! < x,.2 • ■ • < Xi j,_^ and 

v) 0 > x,i - X.2 ^ x,.2 - > . . . > x,.^_j - x,.^ and 
x,.j >0,...x^^ >0; 

vi) 0 < X;! - x.^ < X;2 - X.J, <...< Xj.^_i - X; and 
X;i >0,...Xi;^ >0; and 

vii) 0 > X., - x.2 > X;2 - x,.3 > . . . > x,.^ - x,^^, and 

e) Qa: and o^ju are weighting factors; and 

f) R(x) is a regularizing term. 

The apparatus claimed in claim 21, wherein said processor circuit is 
programmed to produce a graphical representation of said object, for 
controlling a display device. 
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26. The apparatus claimed in claim 25, wherein said processor circuit is 
programmed to produce a time varying graphical representation of said 
object to represent a change of said property over time. 

27. The apparatus claimed in claim 26, further comprising a transmitter for 
5 transmitting said graphical representation to a display device for 

displaying a representation of said property. 

28. The apparatus as claimed in claim 21, further comprising a transmitter 
for transmitting said values representing said property to a display 
device for displaying a representation of said property. 

10 29. The apparatus claimed in claim 28, further comprising a display for 

displaying said representation of said property. 

30. The apparatus claimed in claim 29, wherein said processor circuit is 
programmed to produce successive images representing said property 
at successive instants in time, in response to said representation of 

15 said property. 

31. The apparatus claimed in claim 30, wherein said processor circuit is 
programmed to produce said successive images to depict a three 
dimensional representation of said property in said object which varies 
according to changes in said property over time. 

20 32. The apparatus claimed in claim 21, further comprising a device for 

producing said sets of values representing measurements of said 
property across said object at respective measurement times. 

33. The apparatus claimed in claim 32, wherein said device comprises a 
radioactivity measurement device for measuring radioactivity across 

25 said object. 

34. The apparatus claimed in claim 33, wherein said device includes a 
tomography scanner. 
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The apparatus as claimed in claim 33, wherein said radioactivity 
measurement device comprises a Single Photon Emission Computed 
Tomography (SPECT) imaging device. 

An apparatus for producing a representation of a measurable property 
which varies in time and space, the apparatus comprising: 

a) means for receiving a plurality of sets of values representing 
measurements of said property across said object at respective 
measurement times, each set being associated with a 
respective measurement time; and 

b) means for producing a plurality of sets of values representing 
said property at a plurality of locations throughout said object at 
said respective measurement times, by minimizing a figure of 
merit function relating said values representing measurements 
with said values representing said property at a plurality of 
locations, with a shape constraint imposed on said values 
representing said property at a plurality of locations. 

A computer readable medium for providing computer readable 
instructions for directing a programmable device to implement a 
method of producing a representation of a measurable property which 
varies in time and space by: 

a) receiving a plurality of sets of values representing 
measurements of said property across said object at respective 
measurement times, each set being associated with a 
respective measurement time; and 

b) producing a plurality of sets of values representing said property 
at a plurality of locations throughout said object at said 
respective measurement times, by minimizing a figure of merit 
function relating said values representing measurements with 
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said values representing said property at a plurality of locations, 
with a shape constraint imposed on said values representing 
said property at a plurality of locations. 

A computer data signal embodied in a carrier wave for directing a 
programmable device to implement a method of producing a 
representation of a measurable property which varies in time and 
space, said data signal comprising: 

a) a code segment for directing a programmable device to receive 
a plurality of sets of values representing measurements of said 
property across said object at respective measurement times, 
each set being associated with a respective measurement time; 
and 

b) a code segment for directing a programmable device to produce 
a plurality of sets of values representing said property at a 
plurality of locations throughout said object at said respective 
measurement times, by minimizing a figure of merit function 
relating said values representing measurements with said 
values representing said property at a plurality of locations, with 
a shape constraint imposed on said values representing said 
property at a plurality of locations. 

A method of producing an image representing changes of radioactivity 
in an object, the method comprising; 

a) receiving a plurality of sets of values representing tomography 
scan images across said object at respective measurement 
times, each set being associated with a respective 
measurement time; 

b) producing a plurality of sets of values representing radioactivity 
at a plurality of locations throughout said object at said 
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respective measurement times, by minimizing a figure of merit 
function relating said values representing said tomograpliy scan 
images with said values representing radioactivity, with a shape 
constraint imposed on said values representing radioactivity; 
and 

c) producing a visual representation of said object in response to 
said plurality of sets of values representing radioactivity, said 
visual representation including a representation of radioactivity 
over time. 

An apparatus for producing an image representing changes of 
radioactivity in an object, the apparatus comprising; 

a) a receiver for receiving a plurality of sets of values representing 
tomography scan images across said object at respective 
measurement times, each set being associated with a 
respective measurement time; 

b) a processor circuit programmed to produce a plurality of sets of 
values representing radioactivity at a plurality of locations 
throughout said object at said respective measurement times, by 
minimizing a figure of merit function relating said values 
representing said tomography scan images with said values 
representing radioactivity, with a shape constraint imposed on 
said values representing radioactivity; and programmed to 
produce a time varying graphical representation of said object in 
response to said plurality of sets of values representing 
radioactivity; and 

c) a display responsive to said graphical representation for 
producing a visual representation of said object in response to 
said graphical representation. 



